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Abstract 

The hrst observation of the decays B s — > K~7r + 7r~ and B° — > D~^K~7r + ir~ are 
reported using an integrated luminosity of 1.0 fb -1 recorded by the LHCb experi- 
ment. The branching fractions, normalized with respect to — > D+7r~7r + 7r~ and 
B® — > Df K~7T + 7r~ , respectively, are measured to be 

e ^^ D ' +A '^»=(5.2 ± 0.5 ± 0.3) X 10- 



■> DJ 7T 7T+7T ) 

l _ n i ; = 0.54 ± 0.07 ± 0.07, 

where the first uncertainty is statistical and the second is systematic. The B® — > 
K~ir + 7r~ decay is of particular interest as it can be used to measure the weak 
phase 7. First observation of the B® — > D s i(2536) + 7r~, — > D+tt~tt + decay is 
also presented, and its branching fraction relative to B® — > L>+7r~7r + 7r _ is found to 
be 

IHg-.ft,gyr,mi-.»-,') _ ± l0 ± x 10 _ 3 

■> 7T 7T+7T ) 
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1 Introduction 



In the Standard Model (SM), the amplitudes associated with flavor- changing processes 
depend on four Cabibbo-Kobayashi-Maskawa (CKM) [HI2] matrix parameters. Contribu- 
tions from physics beyond the Standard Model (BSM) add coherently to these amplitudes, 
leading to potential deviations in rates and CP-violating asymmetries when compared to 
the SM contributions alone. Since the SM does not predict the CKM parameters, it is 
important to make precise measurements of their values in processes that are expected 
to be insensitive to BSM contributions. Their values then provide a benchmark to which 
BSM-sensitive measurements can be compared. 

The least well-determined of the CKM parameters is the weak phase 



~ 10° — 12° [21 11]. It may be probed using time- independent rates of decays such as 
B~ — > DK~ [5l El [7j , or by analyzing the time- dependent decay rates of processes such as 
B° s — > DJK ± [HI [91 QUI Ej. Sensitivity to the weak phase 7 results from the interference 
between b — > c and b — > u transitions, as indicated in Figs. [T^a-c). Such measurements 
may be extended to multibody decay modes, such as B~ — >■ DK~tt + ti~ [12] for a time- 
independent measurement, or B° s — > D+K~7t + n~ in the case of a time-dependent analysis. 

The B° — > D + K~7t + tt~ decay, while having the same final state as B Q S — > DfK~7T + 7T~ , 
receives contributions not only from the V^-exchange process (Fig. GJd)), but also from 
b — > c transitions in association with the production of an extra ss pair (Figs. (He-f)). 
The decay may also proceed through mixing followed by a b — > u, ^-exchange process 
(not shown). However, this amplitude is Cabibbo-, helicity- and color-suppressed, and is 
therefore negligible compared to the b — > c amplitude. 

This paper reports the first observation of B° s — > D^K~7t + 7t~ and B° — > D+K~7r + 7T~ , 
and measurements of their branching fractions relative to B° s — > D^7T^7i + ii^ and 
B° s — y DfK~7T + ir~ , respectively. The data sample is based on an integrated luminos- 
ity of 1.0 fb" 1 of pp collisions at ^/i = 7 TeV collected by the LHCb experiment. The 
same data sample is also used to observe the B® — > _D s i(2536) + 7r~ , Df x — > Dfn^7i + decay 
for the first time, and measure its branching fraction relative to B® — > D+tt^tt+tt^ . The 
inclusion of charge-conjugated modes is implied throughout this paper. 

2 Detector and simulation 

The LHCb detector [13] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift-tubes placed downstream. The combined tracking system has 
a momentum resolution (Ap/p) that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and an impact parameter (IP) resolution of 20 |xm for tracks with high transverse mo- 




which, through direct measurements, is known to a precision of 
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Figure 1: Diagrams contributing to the B®,B° S — > D^K~tt + tt~ (a-c) and B° — > 
DfK~7T + 7T~ (d-f) decays, as described in the text. In (a-d), the additional (7r + 7r~) 
indicates that the K~n + ir~ may be produced either through an excited strange kaon 
resonance decay, or through fragmentation. 

mentum (pr)- Charged hadrons are identified using two ring-imaging Cherenkov (RICH) 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. 

The trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage, which applies a full event reconstruction. 
The software trigger requires a two-, three- or four-track secondary vertex with a high 
pr sum of the tracks and a significant displacement from the primary pp interaction 
vertices (PVs). At least one track should have pr > 1.7GeV/c, an IP x 2 greater than 
16 with respect to all PVs, and a track fit x 2 / n df < 2, where ndf is the number of 
degrees of freedom. The IP x 2 is defined as the difference between the x 2 °f the PV 
reconstructed with and without the considered particle. A multivariate algorithm is used 
for the identification of secondary vertices [13] . 
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For the simulation, pp collisions are generated using Pythia 6.4 [T5] with a specific 
LHCb configuration [16]. Decays of hadronic particles are described by EvtGen [17] 
in which final state radiation is generated using Photos [18]. The interaction of the 
generated particles with the detector and its response are implemented using the Geant4 
toolkit [TO] as described in Ref. [20] . 

3 Signal selection 

Signal .B9s decay candidates are formed by pairing a — > K + K~ir + candidate with 
either a 7r~7r + 7r~ (hereafter referred to as X^) or a K~ir + 7r~ combination (hereafter re- 
ferred to as X s ). Tracks used to form the Df and X^ s are required to be identified as 
either a pion or a kaon using information from the RICH detectors, have p^ in excess of 
100 MeV/c, and be significantly detached from any reconstructed PV in the event. 

Signal .D+ candidates are required to have good vertex fit quality, be significantly 
displaced from the nearest PV, and have invariant mass, M(K + K~ir + ), within 20 MeV/c 2 
of the mass [21J. To suppress combinatorial and charmless backgrounds, only those 
candidates that are consistent with decaying through either the <p {M(K + K~) < 
1040 MeV/c 2 ) or K*° (\M(K~ix + ) — m K *o\ < 75 MeV/c 2 ) resonances are used (here, m K *o 
is the K*° mass [21]). The remaining charmless background yields are determined using 
the .D+ mass sidebands. For about 20% of candidates, when the K + is assumed to be a 7r + , 
the corresponding K~ir + 7r + invariant mass is consistent with the D + mass. To suppress 
cross-feed from B° — > D + X decays, a tighter particle identification (PID) requirement 
is applied to the K + in the — > K + K~ir + candidates when |M(il"~7r + 7r + ) — m^+l < 
20 MeV/c 2 (m£)+ is the D + mass J2TJ). Similarly, if the invariant mass of the particles 
forming the candidate, after replacing the K + mass with the proton mass, falls within 
15 MeV/c 2 of the mass, tighter PID selection is applied. The sizes of these mass 
windows are about 2.5 times the invariant mass resolution, and are sufficient to render 
these cross-feed backgrounds negligible. 

Candidate Xj and X s are formed from 7r~7r + 7r~ or K~tt + tt~ combinations, where all 
invariant mass values up to 3 GeV/c 2 are accepted. To reduce the level of combinatorial 
background, we demand that the X^ s vertex is displaced from the nearest PV by more 
than 100 \mv in the direction transverse to the beam axis and that at least two of the 
daughter tracks have p-r > 300 MeV/c. Backgrounds to the — > D^K~n + ix~ search 

from B° s — > Di*^ + 7T~7T + 7r~ or B° s — > DfK~ K + ir~ decays are suppressed by applying more 
stringent PID requirements to the K~ and n + in X s . The PID requirements have an 
efficiency of about 65% for selecting X s , while rejecting about 97% of the favored three- 
pion background. To suppress peaking backgrounds from B° s — > D^D~ decays, where 
.D+ — > 7r + 7r~7r + , A' + 7r~7T + , it is required that M{X^ S ) is more than 20 MeV/c 2 away 
from the mass. 

Signal B meson candidates are then formed by combining a with either an Xj 
or X s . The reconstructed B candidate is required to be well separated from the nearest 
PV with a decay time larger than 0.2 ps and have a good quality vertex fit. To suppress 



3 



remaining charmless backgrounds, which appear primarily in B° — > DfK~7T + 7T~ , the 
vertex separation (VS) x 2 between the D+ and B decay vertices is required to be greater 
than 9. Candidates passing all selection requirements are refit with both mass and 
vertex constraints to improve the mass resolution [22] . 

To further suppress combinatorial background, a boosted decision tree (BDT) se- 
lection [23] with the AdaBoost algorithm [2"4"j is employed. The BDT is trained using 
simulated B° s -»■ D+K~ ~7t + tt decays for the signal distributions, and the high B mass 
sideband in data are used to model the backgrounds. The following thirteen variables are 
used: 

• B candidate: IP %' \ VS x 2 ■, vertex fit x 2 , and pr; 

• Df candidate: Flight distance significance from B vertex; 

• X^ s candidate: IP x 2 ; maximum of the distances of closest approach between any 
pair of tracks in the decay; 

• Xd )8 daughters: min(IP % 2 ), max(IP % 2 ), min(px); 

• Df daughters: min(IP % 2 ), max(IP x 2 ), min(p^), 

where min and max denote the minimum and maximum of the indicated values amongst 
the daughter particles. The flight distance significance is the separation between the 
and B vertices, normalized by the uncertainty. The training produces a single variable, x, 
that provides discrimination between signal decays and background contributions. The 
cut value is chosen by optimizing S(x cut )/ ^/S(x cut ) + B(x cut ), where S(x cut ) and B(x cut ) 
are the expected signal and background yields, respectively, after requiring x > x cut . At 
the optimal point, a signal efficiency of ~90% is expected while rejecting about 85% of the 
combinatorial background (after the previously discussed selections are applied). After 
all selections, about 3% of events have more than one signal candidate in both data and 
simulation. All candidates are kept for further analysis. 

4 Fits to data 

The B° s Df7T~7T + 7T~ and B9s — > DfK~-K + ir~ invariant mass spectra are each modeled 
by the sum of a signal and several background components. The signal shapes are obtained 
from simulation, and are each described by the sum of a Crystal Ball (CB) [23] shape 
and a Gaussian function. The CB shape parameter that describes the tail toward low 
mass is fixed based on simulated decays. A common, freely varying scale factor multiplies 
the width parameters in the CB and Gaussian functions to account for slightly larger 
resolution in data than in simulation. For the B9 — > D^K~n + it~ mass fit, the difference 

between the mean B° s and B° masses is fixed to 87.35 MeV/c 2 |21j . 

Several non-signal 6-hadron decays produce broad peaking structures in the 
D^7t~7t + 7t~ and D+K~h + it~ invariant mass spectra. For B° s — > DfiT~7T + 7T~ , the only 
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significant source of peaking background is from B° s — > D* + tt~7t + tt~ , where the photon 
or 7T° from the D* + decay is not included in the reconstructed decay. Since the full decay 
amplitude for — > D* + 7t~tt + 7t~ is not known, the simulation may not adequately model 
the decay. Simulation is therefore used to provide an estimate for the shape, but the 
parameters are allowed to vary within one standard deviation about the fitted values. 
For £?9% — > DfK~7T + ir~, backgrounds from B^ — > D* + K~7r + 7r~ and from misiden- 

tified B® — > D+7r _ 7r + 7r~ and B° s — )■ D* s + it~ it + it~ decays are considered. The B^ — > 

D* + K~-k + 7t~ shape is fixed to be the same as that obtained for the B° s — > D* + tt~7t + tt~ 
component in the B° s — > D+ir~7T + 7T~ mass fit. This same shape is assumed for both 
B° and where for the former, a shift by the B° — B® mass difference is included. 
For the B Q S — > _D+7r~7r + 7r~ and B Q S — > D* + n^n + n^ cross-feed, simulated decays and 
kaon misidentification rates taken from D* + calibration data are used to obtain their ex- 
pected yields and invariant mass shapes. The cross-feed contribution is about 3% of the 
B° s — > ix~ ix + n~ and B° s — > D* + 7t^7t + tt~ yields; the corresponding cross-feed yields are 
fixed in the B9 g > — > D^K~-K +r K~ fit. The shape is obtained by parameterizing the invari- 
ant mass spectrum obtained from the simulation after replacing the appropriate 7r~ mass 
in Xd with the kaon mass. The combinatorial background is described by an exponential 
function whose slope is allowed to vary independently for both mass fits. 

Figure [2] shows the invariant mass distribution for B° s — > D+n^n+ix^ candidates 
passing all selection criteria. The fitted number of B° s — > D^tt^tt+tt^ signal events is 
5683 ± 83. While it is expected that most of the low mass background emanates from 
B° s — > D* + ir~7T + 7T~ decays, contributions from other sources such as B° s — > D+tt^tt+tt^tt 
are also possibly absorbed into this background component. Figure [3] shows the invariant 
mass distribution for B9s — > D+K~~7t + 7t~ candidates. The fitted signal yields are 402 ±33 

B° D+K~ii + ii~ and 216 ± 21 B° s D+K'tt+tt- events. 

The .D+ mass sidebands, defined to be from 35 to 55 MeV/c 2 on either side of the 
nominal mass, are used to estimate the residual charmless background that may 
contribute to the observed signals. The numbers of B® decays in the sidebands are 
61 ± 16 , Otjj, and 9 ± 5 for the B° s D+tt-tt+tt-, B° s D+R-tt+tt' and B° 
D+K~7t + tt~ decays, respectively; they are subtracted from the observed signal yields 
to obtain the corrected number of signal decays. The yields in the signal and sideband 
regions are summarized in Table [TJ 

5 Mass distributions of X^, s and two-body masses 

In order to investigate the properties of these B^ decays, s Weights [26] obtained from 
the mass fits are used to determine the underlying <, invariant mass spectra as well as 
the two-body invariant masses amongst the three daughter particles. Figure H] shows 
(a) the 7r~7r + 7r~ mass, (b) the smaller n + ir~ mass and (c) the larger n + ir~ mass in 
B° s — > D^7r~7T + 7T~ data and simulated decays. A prominent peak, consistent with the 
ai(1260) _ — > 7r _ 7r + 7r~ is observed, along with structures consistent with the p° in the 
two-body masses. There appears to be an offset in the peak position of the ai(1260) _ 
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Figure 2: Invariant mass distribution for B® — > D^7i^ir + 7i~ candidates. The fitted signal 
probability disrtibution function (PDF) is indicated by the dashed line and the back- 
ground shapes are shown as shaded regions, as described in the text. 

between data and simulation. Since the mean and width of the ai(1260)~ resonance are 
not well known, and their values may even be process dependent, this level of agreement is 
reasonable. A number of other spectra have been compared between data and simulation, 
such as the pp spectra of the Df, X d and the daughter particles, and excellent agreement 
is found. 

Figure [5] shows the corresponding distributions for the B Q S — > D^K~7t + tt~ decay. A 
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Figure 3: Invariant mass distribution for B9n — > D^K~tt + 7T~ candidates. The fitted 
signal (dashed lines) and background shapes (shaded/hatched regions) are shown, as 
described in the text. 



peaked structure at low K~7r + n~ mass, consistent with contributions from the lower-lying 
excited strange mesons, such as the -Ki(1270)~ and i^i(1400)~, is observed. As many 
of these states decay through K*° and p° mesons, significant contributions from these 
resonances are observed in the K~n + and tt + h~ invariant mass spectra, respectively. The 
simulation provides a reasonable description of the distributions in the data. 

Figure E] shows the same distributions for B° — > D+K~7t + tt~ . The K~tt + tt~ invariant 
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Table 1: Summary of event yields from data in the signal and sidebands regions, and 
the background corrected yield. The signal and sideband regions require candidates 
to have invariant mass \M(K + K~-k + ) - m D +\ < 20 MeV/c 2 and 35 < \M(K + K~ti + ) - 
m D + \ < 55 MeV/c 2 , respectively, where rn D + is the mass [2Tj . 



Decay 


Signal Region 


Sideband Region 


Corrected Yield 


B° s ->■ D+TV-n + TV- 
B° s -> D+K-tt + it- 
B° -> D+K-ir+ir- 


5683 ± 83 
216 ±21 
402 ± 33 


61 ± 16 
9±5 


5622 ± 85 
216 ±22 
393 ± 33 




Figure 4: Invariant mass distributions for (a) Xj, (b) smaller ir + 7r~ mass in X^ and (c) 
the larger 7r + 7r~ mass in Xj, from B° s — > n~ n + ix~ decays using sWeights. The points 
are the data and the solid line is the simulation. The simulated distribution is normalized 
to have the same yield as the data. 

mass is quite broad, with little indication of any narrow structures. There are indications 
of K*° and p° contributions in the K~n + and ti~tt + invariant mass spectra, respectively, 
but the contribution from resonances such as the Ki(1270)~ or /^(MOO) - appear to be 
small or absent. In the K~n + invariant mass spectrum, there may be an indication of a 
i^o(1430)° contribution. The simulation, which models the A'~7r + 7r~ final state as 10% 
#i(1270)- 10% ATi(1400)-, 40% K*°ir- and 40% K~p°, provides a reasonable description 
of the data, which suggests that processes such as those in Figs. QJe-f) constitute a large 
portion of the total width for this decay. 

6 First observation of B° 8 -> D sl (2536) + ir- 

A search for excited states, such as D*j —> D^tt~tt + , contributing to the B° s — > 
D+7r~7T + 7T~ final state is performed. Signal candidates within ±40 MeV/c 2 of the nominal 
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Figure 5: Invariant mass distributions for (a) X s , (b) 7r + 7r~ in X s and (c) the K~ir + in 
X s , from B° s — > D+K~ti + ii~ data using sW eights. The points are data and the solid line 
is the simulation. The simulated distribution is normalized to have the same yield as the 
data. 




1000 2000 3000 1000 2000 3000 1000 2000 3000 

K Vit" Mass [MeV/c 2 ] ifjt* Mass [MeV/c 2 ] KV Mass [MeV/c 2 ] 



Figure 6: Invariant mass distributions for (a) X s , (b) tt + 7t~ in X s and (c) the K~ir + in 
X s , from B° — > D+K~7r + iT~ data using sW eights. The points are data and the solid line 
is the simulation. The simulated distribution is normalized to have the same yield as the 
data. 



B® mass are selected, and from them the invariant mass difference, AM = M (D+ n~ n + ) — 
M(D^) is formed, where both 7r + 7r~ combinations are included. The AM distribution 
for candidates in the B° s signal window is shown in Fig. [71 A peak corresponding to the 
-D s i(2536) + is observed, whereas no significant structures are observed in the upper B Q S 
mass sideband (5450—5590 MeV/c 2 ). The distribution is fitted to the sum of a signal 
Breit-Wigner shape convolved with a Gaussian resolution function, and a second order 
polynomial to describe the background contribution. The Breit-Wigner width is set to 
0.92 MeV/c 2 [21], and the Gaussian resolution is fixed to 3.8 MeV/c 2 based on simulation. 
A signal yield of 20.0 ±5.1 signal events is observed at a mass difference of 565.1 ± 
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1.0 MeV/c 2 , which is consistent with the known D sl (2536) + — mass difference of 
566.63 ± 0.35 MeV/c 2 [21]. The significance of the signal is 5.9, obtained by fitting the 
invariant mass distribution with the mean mass difference fixed to 566.63 MeV/c 2 (21], 
and computing — 21n(£ / £ max ) . Here, £ max and C are the fit likelihoods with the 
signal yields left free and fixed to zero, respectively. Several variations in the background 
shape were investigated, and in all cases the signal significance exceeded 5.5. This decay 
is therefore observed for the first time. To obtain the yield in the normalization mode 
(B° s — > Df7T~7T + 7T~), the signal function is integrated from 40 MeV/c 2 below to 40 MeV/c 2 
above the nominal B° s mass. A yield of 5505 ± 85 events is found in this restricted mass 
interval. 



7 Selection efficiencies 

The ratios of branching fractions can be written as 

B(B° S D+K--K+7T-) Y(B° S D+K-tt+tt- 



B(B° S -> D+tt-tt+tt-) Y(B° S D+Tr-ir+Tr- 

and 

B(B° D+K-tt+t:-) Y(B° D+K-n+n- 



X <el (!) 



x 4, x f s /f d , (2) 



B(B° S D+K-n+n-) Y(B° S D+K-tt+tt 

where Y are the measured yields, e* el = e(£>° — > D+n~ix + n~) / e(B® — > K~ix + n~) and 
e rci = e (.B° s ~~ >" DfK^n + 7Y^))/e(B° DfK^n + n~) are the relative selection efficiencies 
(including trigger), and f s /fd = 0.267±0.021 (27] is the B° s fragmentation fraction relative 
to B°. The ratios of selection efficiencies are obtained from simulation, except for the PID 
requirements, which are obtained from a dedicated D* + calibration sample, weighted to 
match the momentum spectrum of the particles that form Xd and X s . The selection 
efficiencies for each decay are given in Table [2J The efficiency of the B° s — > D+7t~7t + tt~ 
decay is about 35% larger than the values obtained in either the B° s — > DfK~7T + 7T~ or 
B° — > D+K~ir + 7T~ decay; the efficiencies of the latter two are consistent with each other. 
The lower efficiency is due almost entirely to the tighter PID requirements on the K~ 
and 7r + in X s . Two additional multiplicative correction factors, also shown in Table [21 
are applied to the measured ratio of branching fractions in Eqs. [T] and [2J The first is a 
correction for the mass veto on M(X d ^ s ), and the second is due to the requirement that 
M(X s> d) < 3 GeV/c 2 . The former, which represents a small correction, is estimated from 
the sWeight-ed distributions of M(Xd, s ) shown previously. For the latter, the fraction of 
events with M(Xd yS ) > 3 GeV/c 2 is obtained from simulation, and scaled by the ratio of 
yields in data relative to simulation for the mass region 2.6 < M(X S ^) < 3.0 GeV/c 2 . 
A 50% uncertainty is assigned to the estimated correction. Based on the qualitative 
agreement between data and simulation in the M(X d)S ) distributions (see Sect. [5]), and 
the fact that the phase space approaches zero as M(Xd yS ) ~ > 3.5 GeV/c 2 , this uncertainty 
is conservative. The relative efficiency between B® — > -D s i(2536) + 7r~, Df ± — > D+7t~tt + 
and B° s D^tt~7t^tt~ is estimated from simulation, and is found to be 0.90 ± 0.05. 
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Figure 7: Distribution of the difference in invariant mass, M{D^tx^ti + ) — M(D+), using 
B° s — > Df7T~7T + 7T~ candidates within 40 MeV/c 2 of the known B° s mass (points) and in 
the upper B® mass sidebands (filled histogram). The fit to the distribution is shown, as 
described in the text. 



8 Systematic uncertainties 

Several uncertainties contribute to the ratio of branching fractions. The sources and their 
values are listed in Table [3j The largest uncertainty, which applies only to the ratio 
!!!fLi^+^ 71 \ n I , is from the b hadronization fraction, f s l f d = 0.267 ± 0.021 127], which 
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Table 2: Selection efficiencies and correction factors for decay modes under study. The 
uncertainties on the selection efficiencies are statistical only, whereas the correction factors 
show the total uncertainty 



Quantity 


B° s -»■ D+Tr-ir+ir- 




B° D+K-n+n- 


Total e (10- 4 ) 


4.97 ±0.08 


3.67 ±0.10 


3.59 ±0.10 


Df veto corr. 


1.013 ±0.003 


1.013 ±0.003 


1.017 ±0.005 


M > 3 GeV/c 2 corr. 


1.02 ±0.01 


1.04 ±0.02 


1.14 ±0.07 



is 7.9%. Another large uncertainty results from the required correction factor to account 
for the signal with M(X Si d) > 3 GeV/c 2 . Those corrections are described in Sect. [7j 

The selection efficiency depends slightly on the modeling of the Xd, s decay. The 
momentum spectra of the B, Df, X^ s and the X^ s daughters have been compared to 
simulation, and excellent agreement is found. The selection efficiency is consistent with 
being flat as a function of M(Xd tS ) at the level of two standard deviations or less. To 
assess a potential systematic uncertainty due to a possible M(A ds )-dependent efficiency, 
the relative differences between the nominal selection efficiencies and the ones obtained by 
reweighting the measured efficiencies by the Xd, s mass spectra in data are computed. The 
relative deviations of 0.5%, 1.1% and 1.2% for B° s -»■ D+A'-tt+tt", B° s -»■ D+tt-tt+tt- 
and B° — > K~ir + n~ , respectively, are the assigned uncertainties. The systematic 
uncertainty on the BDT efficiency is determined by fitting the B° s — > D+tt~7t + tt~ mass 
distribution in data with and without the BDT requirement. The efficiency is found 
to agree with simulation to better than the 1% uncertainty assigned to this source. In 
total, the simulated efficiencies have uncertainties of 1.6% and 1.9% in the two ratios 
of branching fractions. The PID efficiency uncertainty is dominated by the usage of the 
D* + calibration sample to determine the efficiencies of a given PID requirement [22]. This 
uncertainty is assessed by comparing the PID efficiencies obtained directly from simulated 
signal decays with the values obtained using a simulated D* + calibration sample that is 
re-weighted to match the kinematics of the signal decay particles. Using this technique, an 
uncertainty of 2% each on the B° s — > DfK~7T + 7T~ and B° — > D+K~7t + tt~ PID efficiencies 
is obtained, which is 100% correlated, and a 1% uncertainty for B® — > D+7t~7t + 7t~ . The 
trigger is fully simulated, and given the identical number of tracks and the well-modeled 
spectra, the associated uncertainty cancels to first order. Based on previous studies [T2j . 
a 2% uncertainty is assigned. 

The uncertainties in the signal yield determinations have contributions from both the 
background and signal modeling. The signal shape uncertainty was estimated by varying 
all the fixed signal shape parameters one at a time by one standard deviation, and adding 
the changes in yield in quadrature (0.5%). A double Gaussian signal shape model was 
also tried, and the difference was negligible. For the combinatorial background, the shape 
was modified from a single exponential to either the sum of two exponentials, or a linear 
function. For B® — > D+tt~7t + tt~ , the difference in yield was 0.4%. For B° s — > DfK~7T + 7T~ , 
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Table 3: Summary of systematic uncertainties (in %) on the measurements of the ratios 
of branching fractions. 



Source 


B(E°- 


>DfK-ir+-ir-) 


B{B°^D+ K--K+-K-) 


B(B°- 


-^Df-K-iT+ir-) 


B(B°-+D+K-w+tt-) 


fs/ fd 






7.9 


M(X s>d ) > 3 GeV/c 2 




2.2 


7.0 


Efficiency 




1.6 


1.9 


PID 




2.2 


0.0 


Trigger 




2.0 


2.0 


Signal yields 




4.0 


6.9 


Simulated sample size 




3.0 


3.0 


Total 




6.4 


13.4 



the maximum change was 4%, and for B° — > D+K tt + it , the maximum shift was 1%. 
In the B ^ — > D^K~7t + tt~ mass fit, the B9 — > D* + K~ir + 7T~ contribution was modeled 

using the shape from the 5° — > D+7T~7T + n~ mass fit. To estimate an uncertainty from 
this assumption, the data were fitted with the shape obtained from B° s — > D* + K~ir + 7T~ 
simulation. A deviation of 5.5% in the fitted B° — > D^K^tt + tt^ yield is found, with 
almost no change in the B® — > D+K~7r + 7r~ yield. The larger sensitivity on the B° yield 
than the B° s yield arises because these background contributions have a rising edge in 
the vicinity of the B° mass peak, which is far enough below the B° s mass peak to have 
negligible impact. These yield uncertainties are added in quadrature to obtain the values 
shown in Table |3j The uncertainties due to the finite simulation sample sizes are 3.0%. 

The major source of systematic uncertainty on the branching fraction for B° s — > 
.D, 5l (2536) + 7r _ , Df x — > Dfir~7i + , is from the relative efficiency (5%), and on the fraction 
of events with M > 3 GeV/c 2 (10%). This 10% uncertainty is conservatively estimated 
by assuming a flat distribution in M(Xd) up to 3 GeV/c 2 , and then a linear decrease to 
zero at the phase space limit of ~3.5 GeV/c 2 . Other systematic uncertainties related to 
the fit model are negligible. Thus in total, a systematic uncertainty of 11% is assigned to 
the ratio B(B° S D s1 (253Q) + tt~ , D+ DfTT~TT + )/B(B° s D+tt-tt+tt-). 

9 Results and summary 

This paper reports the first observation of the S° — > D^K~tt + it~ , B° — > D^K~7t + tt~ 
and B° s — > _D sl (2536) + 7r~, — > D+ix~n + decays. The ratios of branching fractions are 
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measured to be 



B(B° S -»■ D+K-tt+tt- 

B(B° S D + TT-TT + TT- 
£(5° -»■ D+if-TT+TT- 



£(5° D+K-iT+ir- 
and 

£(5° D + TT-TT+TT-) 



(5.2 ±0.5 ±0.3) x 10' 
0.54 ±0.07 ±0.07, 



(4.0 ± 1.0 ±0.4) x 10" 



where the uncertainties are statistical and systematic, respectively. The B® — > 
DfK~7T + 7T~ branching fraction is consistent with expectations from Cabibbo suppres- 
sion. This decay is particularly interesting because it can be used in a time-dependent 
analysis to measure the CKM phase 7. Additional studies indicate that this decay mode, 
with selections optimized for only 5° — > D^K~ir + 'K~ , can contribute about an additional 
35% more signal events relative to the signal yield in B° s — > DfK^ 1 alone. 

The B° -> D+K~n + ir~ branching fraction is about 50% of that for B° s -»■ D+K~7i + n~ . 
Compared to the B° — > D+K~ decay that proceeds only via a VT-exchange diagram, where 
B(B° -)• D+K~)/B(B° S -> D+K~) ~ 0.1 [21], the ratio B(B° -> D+K-n+v-)/B(B°, -)• 
DfK~7r + ir~) is about 5 times larger. A consistent explanation of this larger B° — > 
D+K~7t + 7t~ branching fraction is that only about 1/5 of the rate is from the ^-exchange 
process (Fig. GJd)) and about 4/5 comes from the diagrams shown in Figs. QJe-f). The 
observed M(X S ), M(K~-k + ) and M(ir + ir~) distributions in Fig. |6]also support this ex- 
planation, as evidenced by the qualitative agreement with the simulation. 
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